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Abstract. We present a novel approach to estimating the intensity mapping signal of any 
CO rotational line emitted during the Epoch of Reionization (EoR). Our approach is based 
on large velocity gradient (LVG) modeling, a radiative transfer modeling technique that 
generates the full CO spectral line energy distribution (SLED) for a specified gas kinetic 
temperature, volume density, velocity gradient, molecular abundance, and column density. 
These parameters, which drive the physics of CO transitions and ultimately dictate the shape 
and amplitude of the CO SLED, can be linked to the global properties of the host galaxy, 
mainly the star formation rate (SFR) and the SFR surface density. By further employing 
an empirically derived SFR — M relation for high redshift galaxies, we can express the LVG 
parameters, and thus the specific intensity of any CO rotational transition, as functions of 
the host halo mass M and redshift z. Integrating over the range of halo masses expected 
to host CO-luminous galaxies, we predict a mean CO(l-O) brightness temperature ranging 
from ~ 0.6 /xK at z = 6 to ~ 0.03 /xK at z = 10 with brightness temperature fluctuations of 
~ 0.1 and 0.005 respectively, at /c = 0.1 Mpc“^. In this model, the CO emission 
signal remains strong for higher rotational levels at z = 6, with {Tco) ~ 0.3 and 0.05 
for the CO J = 6—7-5 and CO J = 10—)-9 transitions respectively. Including the effects of 
CO photodissociation in these molecular clouds, especially at low metallicities, results in the 
overall reduction in the amplitude of the CO signal, with the low- and high-J lines weakening 
by 2-20% and 10-45%, respectively, over the redshift range 4 < z < 10. 
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1 Introduction 

Within the last decade, spectral line intensity mapping has been proposed as an additional, 
complementary probe of the large-scale structure (LSS) of star-forming galaxies during the 
epoch of reionization (EoR) [1-3]. Glimpses into this era have been limited to observa¬ 
tions of individual massive galaxies and quasars at high redshifts, provided by the Hubble 
Space Telescope (HST) and ground-based telescopes. While in the future, we anticipate new 
instruments like Atacama Large Millimeter Array (ALMA) and the James Webb Space Tele¬ 
scope (JWST) providing more detailed views of this “cosmic dawn”, they will be restricted 
by relatively small fields of view and an inability to observe galaxies that are simply too 
faint to detect individually. Intensity mapping offers a complementary glimpse of the three- 
dimensional structure of the high-redshift universe by imaging aggregate line emissions from 
thousands of unresolved objects and studying the large-scale fluctuations in the given line 
intensity due to the clustering of unresolved sources. 

Intensity mapping can be performed using many different spectral lines, the 21 cm neu¬ 
tral hydrogen line being among the most common, given its unique insight into the evolution 
of the neutral IGM during the EoR [4-9]. In this paper, we focus on the millimeter to far- 
infrared rotational transitions of carbon monoxide (CO), a molecule that forms primarily in 
star-forming regions and whose intensity maps promise a wealth of information on the spa¬ 
tial distribution of star formation in the universe [10]. While CO intensity fluctuations have 
already been studied, initially as foreground contaminants to cosmic microwave background 
(CMB) measurements [2] and then as probes of LSS [3, 11-14], these studies have been lim¬ 
ited to the lowest order transitions of the molecule, mainly, CO(l-O) and CO(2-l). These 
lines are often considered because they are typically among the brightest and have redshifted 
frequencies that can potentially be observed from the ground. However, with the advent of 
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ALMA^ and its frequency coverage (84 - 950 GHz), fluctuations in the line emission of many 
high-J CO transitions will potentially be available and can be measured and translated into 
a 3D map of the early universe. Having access to multiple CO rotational lines will further 
facilitate redshift identification and mitigate line confusion, making it possible to statisti¬ 
cally isolate the fluctuations from a particular redshift by cross-correlating the emission from 
different sets of lines [3, 15]. 

Attempts in the recent literature to obtain a theoretical estimate of the CO emission 
signal from high redshift galaxies have relied heavily on empirical relations calibrated from 
local observations and are limited almost exclusively to the ^^CO J=l—)-0 transition line. 
To calculate this mean CO brightness, a simple model is often adopted that connects the 
strength of the CO emission to the abundance of the dark matter halos that host CO luminous 
galaxies. [3] construct this model by first approximating a galaxy’s star formation rate (SFR) 
as a linear function of the mass of the galaxy’s host halo. They then further assume a linear 
relationship between the line luminosity and SFR and adopt the Leo to SFR ratio from 
M 82 to calibrate the proportionality constant. [13] embrace a similar approach, adopting the 
SFR — M relation proposed by [3], but using a set of empirical scaling relations between a 
galaxy’s SFR, far-infrared luminosity, and CO(l-O) luminosity that have been measured for 
galaxies at z < 3. Both studies lead to a simple empirical estimate of the CO luminosity 
that is linear in halo mass {Leo oc SFR oc M) and that relies on the extrapolation of low- 
redshift calibrations to higher redshifts, redshifts corresponding to the EoR. [12] arrives at 
the relation between CO(1-0) luminosity and the halo mass via a different route, making use 
of the Millennium numerical simulation results of [16], a study which, although incorporates 
many physical processes to model the CO emission from high-redshift galaxies, still invokes 
low-redshift measurements to calibrate the normalization factor of the CO luminosity for a 
given halo. 

These various models, among others, have led to estimates of the CO power spectra am¬ 
plitude signal that vary over a range spanning two orders of magnitude, illustrating the lack 
of theoretical understanding of the physics of CO transitions in a high-redshift context [10]. 
In [17], this problem is addressed and a computation of CO fluxes is presented within an 
analytic framework that incorporates both global modes of star formation and the physics 
of molecular rotational lines in z > 6 Lyman-break galaxies. Our paper follows the general 
direction taken by [17] and introduces a simpler approach that captures and ties the physics 
of molecular emission lines to high-redshift {z > 4) observations of star-forming galaxies [18] . 
Our approach is based on large velocity gradient (LVG) modeling, a radiative transfer mod¬ 
eling technique that generates the full CO spectral line energy distribution (SLED) for a 
specihed set of physical parameters. Typically, LVG modeling is employed to quantitatively 
analyze an observed set of emission lines and determine the set of parameters that best re¬ 
produce the observed SED [19]. In this paper, we consider applying the LVG methodology in 
the reverse direction: given a halo of mass M with CO-emitting molecular clouds character¬ 
ized by a kinetic temperature Tkin, velocity gradient dv/dr, gas volume density n, molecular 
abundance X) and column density N, we will derive the full CO SED and compute the mean 
surface brightness of any CO rotational line emitted by halos with that mass at any given 
redshift. 

Our paper is organized as follows. In Section 2, we introduce our LVG model for the 
specihe intensity of CO emission and outline the formalism that relates the set of LVG 
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parameters driving the physics of CO transitions in molecular clouds to the global properties 
of the host galaxy, mainly, the SFR. Given the empirically determined high-redshift SFR — 
M relation [18], these parameters, and thus, the specific CO intensity, can ultimately be 
expressed as functions of the host halo mass M and redshift z. In Section 3 we compute 
the spatially averaged CO surface brightness from star-forming galaxies, as well as the power 
spectrum of spatial fluctuations in the CO emission at any given redshift, with a focus on 
redshifts corresponding to the EoR. We conclude in Section 4 with a summary of our results 
and a brief comparison with other related calculations of the CO intensity mapping signal. 
Throughout we consider a ACDM cosmology parametrized by n* = 1, (Tg = 0.8, 6c = 1.69, 
= 0.31, Ha = 0.69, Qi, = 0.05, and h = 0.7, consistent with the latest measurements from 
Planck [20]. 


2 Modeling the CO Emission 


2.1 CO Brightness Temperature 

To calculate the average CO brightness temperature, we follow the formalism presented 
by [13] and consider the specific intensity of a CO line observed at frequency fobs at redshift 
z = 0, 


li^ohs) - 


dz' 


e[l^obs{^ +z')] 


( 2 . 1 ) 


as determined by solving the cosmological radiative transfer equation, where H{z) is the 
Hubble parameter and e [iZobsi^ + z')] is the proper volume emissivity of the given line. Since 
CO is emitted from within halos hosting star-forming galaxies, we take the CO luminosity, 
Lco: to be some function of the halo mass M and redshift, and assume the profile of each 
CO line is a delta function in frequency. 


Leo = L{M, z)6d{i' - vj) 


( 2 . 2 ) 


where vj is the rest frame frequency of the transition of interest. If we then further assume 
that at any given time, a fraction fduty of halos with mass larger than MminfiO actively emit 
CO lines, then for a given halo mass function dn/dM, the volume emissivity is, 

/■“ dn 

e{i',z) = Soil'- + zf fduty dM —{M,z)L{M,z) (2.3) 

The specific intensity of a line with rest frame frequency uj, emitted by gas at redshift 
zj thus simplifies to 

cl r oo j 

JmJco (2.4) 

or, written as the brightness temperature, 

where ks is the Boltzmann constant. 

To determine L(M,zj), the specific luminosity of a given CO line, we employ large 
velocity gradient (LVG) modeling, a method of radiative transfer in which the excitation 
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and opacity of CO lines are determined by the kinetic temperature Tkim velocity gradient 
dv/dr, gas density n, CO-to-H 2 abundance ratio xcOi and the CO column density of the 
emitting source. We adopt the escape probability formalism [21, 22] derived for a spherical 
cloud undergoing uniform collapse where 

1 — 

Pj = - ( 2 . 6 ) 

TJ 

is the probability, for a line optical depth rj, that a photon emitted in the transition J —?■ J—1 
escapes the cloud (see [23] for a more detailed presentation of the LVG formalism.) Assum¬ 
ing that each emitting source consists of a large number of these unresolved homogeneous 
collapsing clouds, the corresponding emergent intensity of an emission line integrated along 
a line of sight can be expressed as 

Ij = ^^Ajxj^j{tj)Nco (2.7) 

where xj is the population fraction in the level, huj is the transition energy, Aj is the 
Einstein radiative coefficient, and Nco is the beam-averaged CO column density. The LVG- 
modeled specihc luminosity of a line emitted by a host halo with disk radius Rd, therefore 
takes the form 


Lco{M,zj) = ^■k‘^R?^Ij^lvg = '^hujRlAjxjl3j{Tj)Nco ■ (2.8) 


We set 


Rd{M, z) 


A jd 

V2md 


^vir 


A jd 

\^md 


X 1.5 


0.m{z) ISvr^ 


- 1/3 


M 

IO8M0 


V 10 



(2.9) 


where Ac = 187r^-|-82(i—39d^, d = VLrn{z) — '\-, and ^m{z) = VLrn{^+z)^/ (nm(l+-z)^+llA) [15]. 
We assume that the specific angular momentum of the material that forms the disk is the same 
as that of the halo, i.e. jd/rnd = 1, and adopt a spin parameter of A 0.05, corresponding 
to an isolated exponential disk [24]. 

Since the excitation state and optical depth of a given line, xj and tj respectively, are 
determined by the set of physical parameters {Tkin, dv/dr, n, xco} that characterize the 
emitting molecular clouds, the task remains to express these parameters as functions of M 
and z, global properties of the host halo. 


2.2 The Star Formation Model 

As will be physically motivated below, the LVG parameters that dictate the shape of CO 
SLEDs in galaxies are well correlated with the galaxy’s global star formation rate surface 
density. The hrst ingredient of our model is therefore a SFR — M relation that connects the 
SER and host halo mass at the high redshifts we are concerned with in this paper. In [18], such 
an empirical relation is derived by mapping the shape of the observed ultraviolet luminosity 
functions (UV LFs) at z ~ 4-8 to that of the halo mass function at the respective redshifts. In 
this abundance-matching method, each dark-matter halo is assumed to host a single galaxy 
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Figure 1. Left panel: The mean SFR — M relation in the high redshift universe, i.e. z > A, derived 
empirically via abundance matching and fitted by the double power law given in eq. ( 2 . 11 ) where 
{oi, 02 , &i, 62 } = {2.4x10“^^, 1.1x10“®, 1.6, 0.6} for fuv = 1 • Right panel: The corresponding 
SFR surface density, T,sfr{M, z), derived by dividing the SFR by the halo disk area, TTRd{M, at 
redshifts z = 6 (red) and z = 10 (blue). For more details and comparison to data, see [18]. 


and the number of galaxies with star formation rates greater than SFR are equated to the 
number of halos with mass greater than M, 

(] f } ^00 

fuv dM—{M,z)= dSFRcfiSFR,z) (2.10) 

Jm dM Jsfr 

where (j){SFR,z) is the observed, dust-corrected UV LF at redshift z and fuv is the star- 
burst duty cycle, i.e. the fraction of halos with galaxies emitting UV luminosity at any given 
time. [18] find that the SFR — M scaling law remains roughly constant across this redshift 
range and thus, an average relation can be obtained and applied to even higher redshifts, 
z > 8 , where it faithfully reproduces the observed z ~ 9 and 10 LFs. This mean scaling law, 
SFRav{M), is fairly well parameterized by a double power law of the form. 


SFRav{M) 


oiM''! , M <Mc 
a2M^2 , M > Me 


( 2 . 11 ) 


with a turnover at a characteristic halo mass ~ 10^^'® Mq. Fitting the average relations in 
the observed SFR range ~ 0.1 - 500 MQ/yr, we obtain {oi, 02 , 61 , 62 } = {2.4x10“^’^, 1.1x10“®, 
1.6, 0.6} for fuv = 1- We find it reasonable to assume a UV duty cycle of unity throughout 
our calculations given that the time between mergers grows shorter than the Hubble time at 
the high redshifts we are considering and the fact that typical hydrodynamical simulations, 
where star-formation is driven not just by mergers, find a star-forming galaxy in effectively 
every halo at these redshifts [25]. 
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The corresponding mean SFR surface density, T,sfr (units Mq yr“^kpc“^), is computed 
by dividing the SFR from eq. (2.11) by the area of the active star-forming halo disk with 
radius given by eq. (2.9). Figure 1 depicts the average SFR — M relation and the resulting 
SFR surface densities at redshifts z = G and 10. 


2.3 Theoretical Models for LVG Parameters 

The next key step in our LVG-motivated approach to predicting the high redshift CO emis¬ 
sion signal is to model the LVG parameters dictating the shape of the CO SLED as functions 
of the global properties of the host halo. As pointed out in [26], quantities, such as the gas 
temperature and density, which characterize the CO-emitting molecular interstellar medium 
(ISM) are well-correlated with the star formation rate surface densities of galaxies. Qual¬ 
itatively, this makes sense since regions of high SFR density typically arise from denser 
gas concentrations and have large UV radiation fields, with increased efficiency for thermal 
coupling between gas and dust. In the following sections, we will outline how each LVG 
parameter can be expressed in terms of the SFR surface density, F,sfr, and thus ultimately 
as a function of just the halo mass and redshift. 


2.3.1 Gas Kinetic Temperature 

To determine the effective kinetic temperature of the CO-emitting molecular gas, we assume 
that the gas and dust in the star-forming disk are thermally well-coupled, i.e. Tgas ~ T^ust- 
The temperature of a dust grain is set by the balance of radiative heating and cooling 
processes taking place in molecular clouds. The rate of heating due to absorption of optical 
or UV radiation and the incident cosmic microwave background (CMB) can be written as 


dE\ 2 

—— = na 

dt Jabs 


Ljo 


diy Qabs,UV {v)Fv{v) + (tsbTcmb 


( 2 . 12 ) 


where asB is the Stefan-Boltzmann constant, a is the grain radius, Qabs,uvi^) is the emissiv- 
ity in the UV-optical regime, which we set equal to unity [27], Tqmb = 2.73(1-1-2:) is the CMB 
temperature and F^{f) is the flux of energy radiated by the central starburst in the disk. 
Assuming optically thick conditions requires this starlight energy to be totally reemitted in 
the infrared; the integral on the right-hand side can thus be expressed as and, 

adopting the conversion between Lif and SFR presented in [28], 


Fir 

^erg s“^ 

eq. (2.12) takes the final form: 

'2.2 X 10^3 

= 7ra" 

abs 


dE\ 2 

—- = na 
dt , 


= 2.2 X 10 


'Fsfr 


,43 


SFR 


-1 


M© yr 


+ U5b(2.73(1 + 2))' 


dvr (M© yr ^kpc ^) 

The rate of cooling of dust grains by infrared emission is given by 


(2.13) 


(2.14) 


= 4:710^ dv Qabs,IR{T^)T^B^{v,Td) (2.15) 

Jo 

where Qabs,iR{i^) — Qabs{b'o){^/^o)^ is the emissivity in the infrared regime, fq is the reference 
frequency, and fJ is the dust emissivity index. Substituting in for QabsjRi^) and FBy{F,Td) 
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(the flux emitted by a black-body), the integral simplifies to 


fdE\ 2 

V dt Jem 


6 (^ fkBV Qabs<E0) ^^^^ ^ ^ 


(2.16) 


At thermal equilibrium, {dE/dt)abs = {dE/dt)em, and the steady-state grain temperature in 
units of kelvin is 


Td{EsFR{M,z),z)= (-(—') 




1 


GOykBj Qa6s(i^o)r(/3 + 4)C(/3 + 4) 


2.3 X 10 


-3 


^SFR 


{M,z) 


Afusb (MQyr ^kpc ^) 


+ (2.73(1+ z))^ j . 

(2.17) 


1/(/3-H4) 


In the analysis below, we take the fiducial value of /? = 1.3, consistent with the mean emissivity 
index derived from the SCUBA Local Universe Galaxy Survey [29], and the standard value 
(5a6s(125 /rm) = 7.5x10“^ [27, 30, 31]. A plot of the kinetic temperature as a function of 
halo mass at different redshifts can be found in the left panel of figure 2. 


2.3.2 Cloud volume density 


Another key parameter in determining the shape of the CO SLED is n, the cloud volume 
density (cm“^) of the dominant collision partner for CO rotational excitation. Since we are 
assuming CO-emitting molecular clouds, n in this case is the cloud volume density of H 2 , 
given by 


_ 3Sc; 
A^mH2'rci 


(2.18) 


where = 1.36 takes into account the helium contribution to the molecular weight (assuming 
cold, neutral gas), mH 2 = 3.34x10“^^ kg, and rd is the cloud radius assuming a uniform 
sphere. At these high redshifts where the ISM is dominated by molecular gas, the surface 
density of H 2 in a given molecular cloud, Scz, can be related to the beam-averaged gas surface 
density in the galactic disk, Tjgas- 

Empirical studies have found that a correlation exists between the surface density of 
molecular gas and the surface density of the SER, a discovery that is consistent with ob¬ 
servations that stars form predominantly in the molecular component of the ISM. The 
Kennicutt-Schmidt (KS) relation [32, 33] formulates this correlation in terms of a power- 
law, (Esfr) cc {Tigas}^, where estimates of the index N range from super-linear [28, 33-36], 
to linear [37-39], to sublinear [40]. Given our paper’s focus on high redshift star-forming 
sources with CO-emitting molecular clouds, we deemed it most appropriate to adopt the KS 
relation presented in [41], 


^SFR = (3.3 ±0.6) X 10 


-4 


^gas 


1 Mq pc 


-2 


1.2±0.1 


Mq yr ^ kpc 


(2.19) 


a relation that was derived from data sets of CO molecular emission in z ~ 1-3 normal 
star-forming galaxies and restricted to the regime where molecular gas dominates the ISM 
at these redshifts, i.e. Egas > 3 Mq pc“^. Since no evidence of any redshift-dependence of 
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Figure 2. Left panel: Gas kinetic temperature assuming the fiducial values (3 = 1.3 and Qabs{^25 
fim) = 7.5x10“^ in eq. (2.17). Right panel: Cloud volume density of molecular hydrogen as defined 
by equations (2.18)-(2.27) plotted at redshifts z = 6 (red) and z = 10 (blue). 


this relation has been found thus far, applying eq. (2.19) at the higher redshifts we consider 
in this paper, z > 4, is a reasonable extrapolation. 

We therefore set Scz equal to as defined by eq. (2.19), except in cases where Tigas 
drops below the threshold surface density for which the cloud is predominantly molecular. 
This “star-formation” threshold, dehned by a molecular gas fraction of fH 2 = 0.5, is derived 
in [42] for a plane-parallel slab (including H 2 dust) as a function of metallicity Z', 


-‘gas,* 


{Z') = 


2m 


( 1.6 In 

'aG(Z') ' 

“hi 

V 

3.2 


( 2 . 20 ) 


where m = 2.34x10“^^ kg is the mean particle mass per hydrogen nucleus, ag{Z') = 1.9 x 
1Q~‘^^Z' cm“^ is the dust-grain Lyman-Werner-photon absorption cross section per hydrogen 
nucleon, and aG is the dimensionless parameter that defines the LW-band optical depth in 
the cloud due to HI dust. 


aG{Z') 


/i+3.1Z'0'365\ 

V ) 


6.78 

1 -b V2.64Z' ■ 


( 2 . 21 ) 


We adopt the fundamental metallicity relation (FMR), a tight relation between the gas- 
phase metallicity Z\ stellar mass M*, and SFR, to ultimately express Z' as a function solely 
of the halo mass and redshift. The FMR was initially observed and formulated in [43] for 
local galaxies in the mass range 9.2 < logM^/M q < 11.4. Since then, the FMR has been 
confirmed to hold for star-forming galaxies at redshifts as high as z ~ 3 [44], and to extend 
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smoothly at lower masses [45], taking the final form, 


12 + log {O/H) 


8.90 + 0.37m — 0.14s — 0.19m^ + 0.12ms — 0.054s^ for /ro .32 > 9.5 
8.93 + 0.51(/xo.32 - 10) for /ro .32 < 9.5 

( 2 . 22 ) 


where = log (M*) — a log {SFR), m = log (M*) — 10, and s = log (SFR). 

To further parametrize the metallicity as a function of the halo mass and redshift, we 
rely on observations and models that support the conclusion that the SFR in galaxies at 
redshifts z > 4 scales nearly linearly with increasing stellar mass and does not vary by more 


than a factor of order 2 [46] . This behavior is consistent with a crude estimation of the stellar 
mass of a galaxy with a star formation rate SFR: 





dtSFR{t) ~ SFR{z) X tniz) 


(2.23) 


where tniz) is the age of the universe at a given redshift z. Since each halo of interest formed 
at some fraction of the age of the universe, the right-hand side should be multiplied by a 
factor < 1. We therefore calibrate the above expression using the SFR-M* best-fit parameters 
presented in [46] for z ~ 4 - 6 and obtain the following relation 


M*(M, z) = (0.28 ± 0.02) SFRav{M) tniz) (2.24) 


where SFRav{M) is the average SFR for a halo of mass M. We assume this relation continues 
to apply at redshifts z > 6 in the following calculations. 

Armed with the parametrization of SFR introduced in §2.2, and thus a metallicity Z 
expressed solely as a function of halo mass and redshift, we can now return to our model for 
the cloud surface density, T,ci- To ensure the molecular state of each individual cloud, i.e. 
7^2 — Td equal to the beam-averaged gas surface density T,gas{M, z) (derived by 

inverting eq. (2.19)) for all halo masses M > M and floor it to the value T,gas,*{M, z, fH 2 = 
0.5) for all M < M, 


T.ci{M,z) 


^gasiM,z) ifM>M 

Z,gas,* (M, z, fH2 = 0.5) otherwise 


(2.25) 


where M is the halo mass at which 'Zgas{z) drops below T,gas,*ifH2 = 0.5, z) at a given 
redshift z. 

The other variable that appears in our definition of the H 2 volume density is r^, the 
molecular cloud radius. Assuming that the molecular gas within the clump is in hydrostatic 
equilibrium, the radius of the cloud can be related to its surface mass density through the 
relation, 


rd = 


^s,eff 


(2.26) 


where G is the gravitational constant and Cg^eff is the effective sound speed in the gas 
(taking into account turbulence, -|- (^turbulence) ^ which we set to 10 km/s. The final 

expression for the molecular cloud volume density then simplifies to 


nH2 {M, z) 


SttG 

4.limH2cl^ff 




(2.27) 
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Figure 3. Left panel: Velocity gradient dv/dr as a function of halo mass M assuming self-gravitating, 
virialized clouds as defined by eq. (2.28). Right panel: CO-to-H 2 abundance ratio via eq. (2.29), 
plotted at redshifts z = 6 (red) and z = 10 (blue). 


a plot of which can be found in the right panel of figure 2. The steep decline in number 
density at low halo masses, as depicted in figure 2, mirrors the steeply declining SFR — M 
relation at the low-mass end (see left panel in figure 1); as the star-formation rate diminishes 
by several orders of magnitude with decreasing halo mass, the SFR surface density drops 
accordingly and ultimately translates into reduced gas column and number densities at these 
low masses. 


2.3.3 Velocity Gradient 


For the sake of simplicity, we assume self-gravitating, virialized molecular clouds, in which 
case, the velocity gradient dv/dr and cloud volume density nH 2 are related in the following 
way [47] 


dv 

— ~ 3.1 
dr 


nH2 


104 


cm 


-3 


kms ^pc ^ 


(2.28) 


where is defined in eq. (2.27). A plot of dv/dr as a function of halo mass at different 
redshifts can be found in the left panel of hgure 3. 


2.3.4 CO-to-H 2 Abundance Ratio 

Studies have shown that at high metallicities, i.e. Z' > 10“^, the dominant metal-bearing 
molecule in the ISM is CO. Furthermore, [48] find that even at low metallicities, 30-100% of 
the available carbon is always locked in CO if the CO is shielded, provided that the hydrogen 
gas is in molecular form. In this limit, the relative CO-to-H 2 abundance varies approximately 
linearly with metallicity [48], and assuming most of the carbon is in fact locked up in CO, 
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Figure 4. Beam-averaged CO column density derived at redshifts z=6 (red) and z=10 (blue) under 
conditions where CO photodissociation is accounted for (solid curves, eq. (2.30)) and neglected 
(dashed curves, eq. (2.32)). In the case where the effects of CO photodissociation are included, the 
CO column density drops to zero for halos smaller than M < 10^° Mq, indicating that the CO in the 
gas has fully disassociated. 


XCO is given by, 

Xco(^) - 3x 10-^Z' (2.29) 

where an expression for the relevant metallicity Z' can be found in eq. (2.22). A plot of xco 
is shown in the right panel of figure 3. 

2.3.5 CO Column Density, including photodissociation 

The CO column density, which sets the overall scaling and amplitude of the CO SLED, can 
be expressed most simply as the product of the CO-to-H 2 abundance ratio and the beam- 
averaged molecular column density, 


Nco = XcoNh2 


(2.30) 


where Nh 2 is a power-law function of the SFR surface density, derived by inverting the KS 
relation in eq. (2.19), 


Nh2{^sfr{M, z)) 


2 X 10-^ / ^SFR 

IJ.mH2 VI M© yr-i kpc-V 


(2.31) 


The above expression for Nco does not account for conditions under which CO has 
photodissociated into C and 0"*“ while the gas continues to remain molecular due to either 
H 2 self-shielding or dust-shielding. Such conditions, which may exist on the surfaces of 
molecular clouds or the clumps contained within such clouds, ultimately result in a fraction 
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of H 2 gas that is “dark” in CO transitions. Observations indicate that the column density 
of this “dark gas” can be as high as 30% (~ 3x10^^ cm“^) of the total molecular column 
density in the local Galaxy cm“^) [49]. A first-order approximation of the CO column 

density that accounts for the effects of CO photodissociation is then 

Nco = Xco ( A'j /2 -(2.32) 

where, as before, Z' is the metallicity in solar units. The second term in the above expression 
is essentially the H 2 column required to enable CO dust-shielding under the assumption that 
the threshold for CO survival is set by a universal dust opacity. This approximation captures 
the effects of decreasing metallicity on the abundance of CO, reducing the column density as 
the number of shielding dust grains grows sparse. In the limiting case where Z' drops so low 
that there is not enough dust to effectively shield CO, (corresponding to the parenthesized 
portion of eq. (2.32) growing negative), the CO is fully dissociated and Nco is set to zero. 
We note that eq. (2.32) is only a first-order approximation and that the constant, 3x10^^, 
can be altered due to variations in UV field intensity or gas clumping factors. The beam- 
averaged CO column density at different redshifts is plotted in figure 4, both for the case 
where CO photodissociation is accounted for (solid curves) and neglected (dashed curves). 

2.4 Model CO SLEDs 

Equipped with analytic expressions for the LVG parameters that dictate the shape and 
magnitude of the CO SLED, we can now compute the intensity of each CO line and the 
resulting SED generated by the molecular clouds in a halo with mass M at redshift 2 . To carry 
out the computations, we use the Mark &: Sternberg LVG radiative transfer code described 
in [50], with GO-H 2 collisional coefficients taken from [51] and energy levels, line frequencies, 
and Einstein A coefficients taken from the Gologne Database for Molecular Spectroscopy 
(CDMS). For a given set of parameters, {T^in, dv/dr, xco-, the code determines 

the level populations by iteratively solving the equations of statistical equilibrium which 
balance radiative absorptions, stimulated emission, spontaneous emission, and collisions with 
H 2 using the escape probability formalism discussed in §2.1. Once the level populations are 
computed, the full GO rotational ladder and line intensities follow from eq. (2.7). 

Figure 5 shows the GO SLEDs generated by halos at redshift z = 10 with masses in the 
range M = 10®-10^^ Mq when the effects of CO photodisassociation are excluded. In the 
case where photodissociation is considered, the line intensities emitted by low-mass halos with 
M < 10^*^ Mq (red and orange curves) entirely disappear, reflecting the full dissociation of 
CO in these molecular clouds where dust-shielding has grown inefficient. However, the other 
curves, corresponding to line emission from higher mass halos that have been normalized to 
the ground state, remain unchanged. This is due to the fact that “turning on” photodisso¬ 
ciation merely reduces the beam-averaged CO column density according to eq. (2.32); since 
Nco controls the overall amplitude of the SLED (and not the shape), adjusting this quantity 
simply amplifies or reduces the intensity of all the CO lines by the same amount, leaving the 
ratio between lines unchanged. 

As expected, we find that as the physical conditions in the emitting molecular clouds 
grow more extreme, the CO rotational levels become increasingly populated and the SLED 
rises accordingly. Consequently, the line intensities not only grow in magnitude, but the peak 
of the CO SLED also shifts to higher J values, reflecting the excitation of the more energetic 
states of the molecule. 
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Figure 5. CO SLEDs generated by a halo at redshift z = 6 (dashed curves) and z = 10 (solid 
curves) with halo mass ranging from 10® to 10^® Mq, neglecting the effects of photodissociation. If 
CO photodissociation is taken into account, the red and orange curves, corresponding to line emission 
from halos with M < 10^° Mq, would disappear while the other curves would remain unchanged. 
The brown curve denotes the expected SLED in the optically-thick case when the population levels 
are in local thermal equilibrium (LTE); in this limiting case, the intensity at a given level follows the 
Planck function, B{v,T)dv = jdz/, where (T) ~ 91 K. 


The curves plotted in figure 5 demonstrate this trend in the case where a UV duty cycle 
of unity is assumed. In this scenario, the SFR and SFR surface density range from 10“^- 
1000 MQyr“^ and 10“^-20 MQyr“^kpc“^ respectively, when the halo mass varies from 10® 
to 10^® Mq (see figure 1). While the corresponding kinetic temperature remains relatively 
constant for this T^sfr range {Tkin ~ 90 K), the H 2 number density varies significantly, 10^ 
cm“® < 71^2 ^ 3 X 10^ cm“®. The variance in the shape of the CO SED computed for 
fuv = 1 reflects this range of physical conditions parameterized by the halo mass at z = 
10; the SEDs produced by low-mass halos, M < 10^*^ Mq, peak at J ~ 4 before turning 
over and plummeting (neglecting CO photodissociation). In these low density clouds, with 
correspondingly small velocity gradients {dv/dr ~ 0.5 kms“^pc“^) and CO abundances 
ixco ~ 5x10“®), only the low-lying transitions such as CO J = 1 —)• 0 are optically thick. In 
contrast, in more extreme star-forming galaxies, i.e. M > 10^® Mq, the H 2 number densities 
reach ~ 10® cm“®, the typical critical density value for high-J CO emission. Consequently, the 
high-J lines grow optically thick and the line ratios approach thermalization for transitions 
as high as J ~ 13 in these high-mass halos. 
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Figure 6. The CO rotational line flux as a function of halo mass at z = 6 (left panel) and z = 10 
(right panel) for Jupper = 1-8. Solid and dashed curves denote results obtained by neglecting (eq. 
(2.30)) and including (eq. (2.32)) the effects of CO photodissociation. 


3 Results 


3.1 Predicted CO Fluxes and Spatially Averaged CO Brightness Temperature 

Given our LVG model of the full GO SLED, we can now predict the GO line fluxes generated 
by a set of molecular clouds residing in a host halo of mass M at redshift z, as well as 
compute the spatially averaged brightness temperature of any GO rotational transition. The 
former is obtained by converting LVG-derived GO luminosities (eq. (2.8)) into observed, 
velocity-integrated fluxes, Fco, using typical observer units. 


= 1.040 X 10“^ 
L© 


Dl \^ T^obs Fco 
Mpc ) GHz Jy km s“^ 


(3.1) 


where is the luminosity distance. The results are shown in figure 6 for a range of GO 
rotational lines emitted as a function of host halo mass at z = 6 (left panel) and z = 10 
(right panel). The dashed curves represent the fluxes obtained by including the effects of GO 
photodissociation in our computations; as expected, the change in flux when accounting for 
this phenomenon is most pronounced at the low-mass end where Fco drops to zero once the 
GO is fully photodissociated. At the higher-mass end, we predict that a ~ 10^^ M© halo at 
redshift z = 6 will emit strongest in the J = 7 —)• 6 transition with a flux of ~ 25 mJy while 
a halo with the same mass at z = 10 will emit strongest in a higher energy state, J = 10 —?• 
9, but with nearly half the flux, ~ 14 mJy. 
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Figure 7. Volume-averaged CO brightness temperature as a function of redshift for a minimum host 
halo mass of CO luminous galaxies of M^m,co = 10 ® (red), 10 ® (green), and 10 ^® Mq (blue curves), 
neglecting the effects of CO photodissociation. The dashed curves denote the signals obtained when 
photodissociation is taken into account {Mmin,co = 10 ^®). Each panel shows {Tco) for a different 
line in the CO rotational ladder with 1 < Jupper < 8. 


In order to determine the spatially-averaged CO brightness temperature emitted by 
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halos across the mass range, we compute the following integral, 


where vj is the rest-frame frequency of the specified line. This equation is parameterized by 
Mmin,co, the minimum host halo mass for CO luminous halos, and fduty, the duty cycle for 
CO activity. Since CO lines are excited by starburst activity, we generally expect that the 
duty cycle for CO luminous activity is comparable to the starburst duty cycle. We therefore 
assume fduty = fuv = 1 in our fiducial models. Furthermore, in computing the volume- 
averaged CO brightness temperature, we vary Mmin,co widely between Mmin,co = 10 ®, 10 ^ 
and 10^^ Mq to illustrate the sensitivity of the results to this parameter. In our model, these 
halos host CO luminous galaxies with SFRs of 10 3x10 and 0.1 M 0 yr ^ (left panel 

of figure 1 ). 

We plot the results in figure 7, where the mean brightness temperature is shown as a 
function of redshift for different CO rotational lines. As expected, we find that (Tco) is a 
steeply declining function of z, a direct consequence of the decreasing number of host halos 
(per volume) at these high redshifts predicted by the Sheth-Tormen halo mass function [52] . 
When CO photodissociation is neglected (solid curves), the redshift evolution of (Tbo) for 
the low-J CO lines steepens at the high-z end when larger values are assumed for Mmin,co- 

While halos with masses M > 10^*^ Mq emit the most flux at all redshifts (figure 6 ), the 
flux emitted by halos with M ~ 10^-10^® Mq makes a non-negligible contribution to the total 
CO brightness temperature at higher redshifts. An illustration of this is presented in figure 8 
which plots dT/dinM, the contribution to the mean brightness temperature per logarithmic 
halo mass for different CO lines as a function of halo mass at 2 : = 6 (dashed curves) and 2 = 
10 (solid curves). It is clear that at z = 10, although the primary component of the CO signal 
originates from halos with masses in the range M lO^'^-lO^^ Mq, the CO emission from 
~ 10® Mq halos makes up nearly 15% of the total signal for the low-J lines. Thus, raising 
the minimum host halo mass for CO luminous halos from 10® Mq to 10^® Mq excludes a 
population of CO-emitting sources and reduces the amplitude of {Too) for the low-J lines by 
a factor of ~ 2. The effects of raising Mmin,co vanish in the higher energy states since these 
low-mass halos emit less than 1% of the total signal in these high-J lines; hence, the red, 
green, and blue solid curves, denoting the mean brightness temperature for models where 
Mmin,co = 10®, 10®, and 10^® Mq, respectively, are barely distinguishable from one another 
for the higher energy rotational transitions. 

In the case where CO photodissociation is taken into account, we found that the CO 
becomes fully photodissociated in halos with M < 10^® Mq (figure 4); since these low-mass 
halos do not emit CO flux in this model, the minimum host halo mass of a CO luminous 
galaxy is effectively set to Mmin,co = 10^® Mq. Accounting for CO destruction in molecular 
clouds results in an overall reduction in the amplitude of (Tco), with the low- and high-J 
lines weakening by ~ 2-20% and 10-45%, respectively, over the range of redshifts shown in 
figure 7 (dashed curves). 

The spatially averaged brightness temperature of the CO J = 1—?- 0 line predicted by 
our model assuming Mmin,co = 10® Mq ranges from ~ 0.6 /uK at z = 6 to ~ 0.03 at z 
= 10 when CO photodissociation is neglected. The strength of the CO signal dwindles for 
higher-J lines, with {Tco{z = 6 )) ~ 0.3 /iK and Tco{z = 10)) ~ 0.1 for the CO J = 6 —)■ 5 
transition. “Turning on” photodissociation does not significantly effect the mean brightness 
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Figure 8. Contribution to the mean brightness temperature per logarithmic mass of different CO 
lines at redshift z = 6 (dashed) and z = 10 (solid), in the case where the effects of CO photodissociation 
are neglected. 


temperature of the high-J lines, but it does reduce the CO(l-O) signal to (Tco) ~ 0-4 and ~ 
0.01 /rK at z = 6 and 10, respectively . 

Without delving into any particular instrumental design, we briefly consider the plausi¬ 
bility of detecting such signals. We use the radiometer equation for the signal-to-noise ratio, 
S/N = [Tco/Tsys)^tint where Tgys ~ 20 K is the system temperature of the detector, 
Az/ is the observed bandwidth, and tint is the integration time, which we take to be 1000 
hours. We find that at ^ = 6, the predicted brightness temperatures of the 4 lowest-lying 
CO transitions (neglecting CO photodissociation) can be detected at 5a confidence with a 
bandwidth of Av ~ 10 GHz, while a bandwidth of Au ~ 100 GHz is required for detection 
with 1(T confidence at z = 10. Higher J lines in this model will require even longer integration 
times to achieve the desired brightness sensitivity. 

3.2 CO Power Spectrum 

Given the existence of brighter foreground sources of emission at the relevant frequencies, the 
mean redshifted CO signal will be difficult, if not impossible, to directly observe. We therefore 
consider spatial fluctuations in the surface brightness and compute the CO power spectrum 
predicted by different variations of our model. In contrast to the spectrally smooth foreground 
sources, the CO signal is expected to have structure in frequency space which can be used to 
isolate spatial fluctuations in its brightness temperature. Since the power spectrum captures 
the underlying matter distribution and structure, a map of the CO brightness temperature 
fluctuations at z > 6 can be used to probe the spatial distribution of star-forming galaxies 
during the EoR. 
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Following the formalism in [12] and [13], the three-dimensional power spectrum of the 
CO brightness temperature fluctuations is expected to take the form 

Pcoik,z) = {Tco){z) [bcoizfPiin{k,z) + Pshotiz)] . (3.3) 


Since CO is emitted from within halos, the first term in this expression represents spatial 
variations due to correlations with the underlying dark matter density field where Pun is the 
linear theory density power spectrum and the bias bco is given by 


bco{z) 


r 

Ji 


M„ 


dM^Lco{M,z)b{M,z) 


'i,CO 


L 


M„ 


dM^Lco{M,z) 


.,co 


(3.4) 


where b{M, z) = l + z) — l)/(ic, I'iM, z) = 5cIcf{M, z), and cr(M, z) is the RMS density 

fluctuation in a spherical region containing mass M [53]. The second term of eq. (3.3), the 
shot noise contribution due to Poisson fluctuations in the number of halos on the sky, can be 
expressed as 

„ , , 1 fS„WALco(M,, 

Pshotiz) = - --- 

iSZ,J¥ml^co(M.z 

where, in all of our calculations, we adopt the Sheth-Tormen halo mass function for dn/dM. 

The power spectra of the CO rotational lines from Jupper = 1 to Jupper = 10 are plotted 
in figure 9 for models which include (dashed) and exclude (solid) the effects of CO photodis¬ 
sociation. The y-axis shows A^q(A:, z) = k^Pcoik, z)/{ 27 r‘^), the contribution to the variance 
of {Too) per logarithmic bin, in units of [/rK^]. In general, the fluctuations depend strongly 
on the wavenumber, with the overall shape of the predicted power spectra reflecting the form 
of eq. (32). On large scales, small k, the clustering term dominates and the fluctuations in 
CO brightness temperature mirror the underlying dark matter density field. Conversly, on 
small scales (large k), the shot noise takes over and the fluctuations reach ~ 300 - 1300 
on scales of /c = 10 Mpc“^ at z = 6, and A^q ~ 4 - 14 /rK^ at z = 10 (when pho¬ 
todissociation is neglected). The J = 1—)-0 transition is predicted to produce the strongest 
signal, with an amplitude that drops by a factor of ~5-50 across the wavenumber range k ~ 
10“^-10 Mpc“^ when the higher energy state, J = 10—>-9, is considered instead. 

The redshift evolution of A^q is ultimately dictated by the behavior of {Tco)iz)', al¬ 
though increases as the host halos become more clustered at higher redshifts, this effect is 
not enough to compensate the declining brightness temperature with z. The auto-correlation 
signal of all the CO lines therefore weakens as the redshift varies from z = 6 to 10, dropping 
by rsj 2 orders of magnitude over this redshift range. On the other hand, our results for A^q 
depend very weakly on our choice of the minimum host halo mass for CO luminous galaxies; 
therefore, in the fiducial models where photodissociation is neglected, we only plot results 
for the case where Mmin,co = 10® Mq, thereby including the contribution of the low-mass 
halos to the overall emission signal. 

The effects of including CO photodissociation are most prominent for the low-J lines, 
weakening the signal by, at most, a factor of ~ 6 at ^; = 10 on large scales. At A: = 0.1 Mpc“^, 
we find CO(l-O) brightness temperature fluctuations of amplitude A^q ~ 0.1 /rK^ at z = 6, 
whether or not CO photodissociation is taken into account. At 2 ; = 10, A^q(A: = 0.1 Mpc“^) 
~ 5x10“^ and 8x10“^ /rK^ for the J = 1—t-O line in the models where photodissociation is 
turned “off’ (assuming fixed Mmin,co = 10® M©) and “on”, respectively. 
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Figure 9. Auto power spectrum of CO brightness temperature fluctuations for lines all the way 
up to Jupper = 10, at z = 6 (red) and z = 10 (blue). In each panel, the solid curves represent 
results for the models in which CO photodissociation was neglected and a minimum host halo mass of 
Mco.min = 10® Mq was assumed; dashed curves denote results for the case where the effects of CO 
photodissociation were included, effectively setting the minimum host halo for CO luminous galaxies 
to Mco.min = 10 ^° M 0 . 


4 Discussion 

We have presented a new approach to estimating the mean CO emission signal from the 
epoch of reionization (EoR) that links the atomic physics of molecular emission lines to 
high-redshift observations of star-forming galaxies. This method is based on LVG modeling, 
a radiative transfer modeling technique that generates the full CO SLED for a specified 
set of characterizing parameters, namely, the kinetic temperature, number density, velocity 
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gradient, CO abundance, and column density of the emitting source. We showed that these 
LVG parameters, which dictate both the shape and amplitude of the CO SLED, can be 
expressed in terms of the emitting galaxy’s global star formation rate, SFR, and the star 
formation rate surface density, Employing the SFR-M relation empirically derived 

for high-redshift galaxies, i.e. z > 4, we can then ultimately express the LVG parameters, 
and thus, the specihc intensity of any CO rotational line, as functions of the host halo mass 
M and redshift z. 

Adopting a starburst duty cycle of fjjv = 1; the average SFR — M relation derived 
via abundance-matching for 4 < z < 8 is characterized by a steeply declining slope at the 
low-mass end, where the star formation rate goes as SFR oc With the SER dropping 

to values below 0.1 MQyr“^ for M < 10^*^ M©, the physical conditions in these low-mass 
halos are not “extreme” enough to substantially excite the high-J CO rotational states. The 
resulting CO SLEDs correspondingly peak around J ~ 4 before turning over, and the overall 
contribution of CO emission from this halo population grows negligible at lower redshifts. 
On the other hand, the H 2 number density and CO-to-H 2 abundance in halos with masses 
M > 10^*^ Mq are large enough to keep the high-J population levels thermalized; the CO 
SLEDs generated by these massive halos therefore have peaks shifited to J > 10 and overall 
higher amplitudes, reflecting the excitation of the more energetic states of the molecule. 

We also consider the effects of CO photodissociation on the CO line intensities generated 
by our fiducial models. Assuming that the threshold for CO survival is set by a universal 
dust opacity, we adopt a first-order approximation of Nco which effectively reduces the CO 
column density with decreasing metallicity. In the limiting case where the metallicity drops 
so low that there is not enough dust to shield CO, the CO is considered fully dissociated and 
Nco is set to zero. We find that such conditions occur in halos with M < 10^*^ Mq, causing 
the line intensities emitted by these low-mass halos to entirely disappear in models where 
CO photodissociation is accounted for. 

Given our LVG model of the full GO SLEDs, we can predict both Fco^ the CO line 
flux generated by a set of molecular clouds in a host halo of mass M at redshift z, as well as 
{Tco)i the spatially averaged brightness temperature of any CO rotational transition. We 
find that the flux emitted in the J = 1 —)• 0 transition by a halo at redshift z = Q with mass 
M 10^^ Mq is Fco{i-o) ~ 0.8 mJy km/s. The higher rotational lines are expected to be 
even brighter, with a 10^^ Mq halo at z = 6 emitting a CO(6-5) flux of ~ 23 mJy km/s. 
These fluxes drop by 25-30% when considering a 10^^ Mq halo emitting at z = 10, with 
^CO(i-o) ~ 0-2 mJy km/s and Fco{e- 5 ) ~ 7 mJy km/s. 

The CO line fluxes emitted by individual host halos of mass M at redshift z as estimated 
in this paper, are found to be generally higher than previous estimates obtained in the 
literature [16, 17]. Building an analytic formalism within a paradigm where star formation 
is a function of gas supply and stellar feedback, the model presented in [17] provides radial 
distributions of SER, 'F,gas, and T which are then used to calculate the masses, sizes, and 
number of GMCs, and the corresponding total CO line luminosities emitted by these clouds. 
Consequently, the CO(l-O) flux emitted by a 10^^ Mq halo at z = 6 as predicted by this model 
is ~ 0.01-0.03 mJy km/s, an order of magnitude smaller than the flux estimates derived with 
our methodology. The flux in the J = 6 —)• 5 transition is found to be highly sensitive to 
the inclusion of turbulent clumps within the GMCs, given that the high densities in such 
inhomogeneities can effect the thermalization of the higher-J rotational lines; the resulting 
CO(6-5) flux in the models presented in [17] thus varies from 10 mJy km/s to 10“^ mJy 
km/s when turbulent clumps are included and excluded, respectively. 
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While these results are comparable to those derived using the semi-analytic methods 
of [16], they are consistently smaller than the values presented in this paper. The prolific 
number of parameters used to characterize the GMC properties in [17] make it difficult to 
ascertain the source of divergence between our results, which both use radiative transfer tech¬ 
niques. However, we suspect these differences can be traced back to the different prescriptions 
used to set the star formation rate surface density, T^sfr, and gas surface density, the 

two ingredients which, once related to one another via the Kennicutt-Schmidt relation in 
our model, determine the shape and amplitude of the resulting CO SLED. Furthermore, [17] 
parameterizes the properties of MCs as a function of their position relative to the disk center, 
and after considering the level populations and optical depth of CO in each cloud directly, 
computes the resulting flux from the entire galaxy by counting up the number of clouds at 
each galactic radius. The molecular clouds in our model, on the other hand, are characterized 
by disk-averaged properties of the host halo disk; to derive the signal from a CO-emitting 
galaxy, we assume a large number of these identical homogeneous collapsing clouds and scale 
the line intensities with the disk-averaged CO column density. Based on the sensitivity lim¬ 
its quoted in [17], our approach predicts that the Jupper = 1 line in z = 6 halos with mass 
M > 10^^ Mq will be observable by JVLA^ (Jansky Very Large Array) after ten hours of 
observation; at this redshift, the CO(2-1) and CO(3-2) lines fluxes emitted by halos with 
M > 10^^ Mq are also expected to be observed by JVLA. Similarly, we expect ALMA to 
detect the CO(6-5) flux emitted hy z = 6 halos with mass M > 5x10^^ Mq after ten hours 
of observation, and higher rotational lines Jupper > 7 emitted by halos with M > 10^^ Mq. 

To obtain an estimate of the spatially averaged brightness temperature of a given line at 
a particular redshift, we simply integrate Lco{M, z) over the range of halo masses that are 
expected to host CO-luminous galaxies. In the case where CO photodissociation is included, 
the minimum host halo mass of CO-emitting galaxies is set to Mco,min ~ 10^*^ Mq by the 
model itself, since the CO is found to be fully dissociated in halos with M < 10^*^ Mq. 
When CO photodissociation is ignored, varying Mco,min from 10® to 10^° Mq reduces the 
low-J line signals at high redshifts, where the CO emission from ~ 10® Mq halos make up a 
non-negligible percentage of the total emission. 

In our fiducial model where CO photodissociation is neglected and Mmin,co = 10® Mq, 
we predict a spatially averaged brightness temperature of {Tco) ~ 0.5 fiK at z = 6 and 0.03 
/iK at z = 10 for the low-J CO rotational lines, with brightness temperature fluctuations of 
amplitude ~ 0.1 and 0.005 /rK^ respectively, at k = 0.1 Mpc“^. These CO emission 
signals are further reduced to {Tco) ~ 0.4 and 0.01 at z = 6 and 10, respectively, for the 
low-lying states when the effects of CO photodissociation are included in the calculations. 
(Note that, since CO lines are typically excited by starburst activity, the choice of Mmin,co = 
10® Mq is favored by theoretical and numerical investigations which indicate that 10® Mq is 
the minimum mass required for a halo to cool and form stars at these high redshifts [15, 54- 
56]). Our estimates of {Tco) for the low-J CO transitions are comparable to the values 
obtained in previous work by [11], [12], [13], and [14]. Constructing a model based on the 
required cosmic star formation rate density to reionize the universe, [11] obtains an order-of- 
magnitude estimate of {Tco){z = 8) ~ 1 /rK for the J = 1—)- 0 and J = 2—)-1 transitions. [12] 
arrives at a slightly smaller estimate of the CO(l-O) brightness temperature, {Tco) ~ 0.5 
/rK for z = 6 and 0.1 /uK for z = 10, by using the Millennium numerical simulation results 
of [16] to model the CO emission from high-redshift galaxies. [13] assumes a linear SFR — M 

^http://www. via.nrao.edu/ 
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relation and a set of \ow-z empirical scaling relations between a galaxy’s SFR, LpiR, and 
Lco(i-o) to estimate {Tco) for these low-J lines; they find a mean brightness temperature 
of ~ 2 at z ~ 6 and 0.5 fiK at z ~ 10 with fluctuations ~ 0.2 and 0.02 
respectively, on scales of A; ~ 0.1 Mpc“^. In a more recent paper, [14] follows the approach 
taken in [13] with a few adjustments to the SFR — M prescription to arrive at brightness 
temperatures of {Tco) ~ 0.7 and 0.2 /iK at redshifts z = 6 and 10 respectively, assuming 
Mmin,co = 10® Mq for the CO(l-O) and CO(2-l) lines. 

While these previous works are limited almost exclusively to predicting the signals of 
the J=l—)-0 and 2—)-l transition lines, our LVG-based approach generates the full CO 

SLED and thus allows us to compute the signal strength of the higher-J energy states as well. 
For example, we predict a CO(10-9) brightness temperature of {Tco) ~ 0.05 at z = 6 and 
0.003 at z = 10 with A‘^Q{k = 0.1) ~ 0.003 and 10“® respectively. We look forward 
to future experiments, such as the Carbon MonOxide Mapping Array (COMA)^ currently 
under development, which promise to provide spectral-spatial intensity mapping of CO at 
the high redshifts characterizing the epoch of reionization. 
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